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Radial velocities of five globular clusters obtained with AAOmega 
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Using the recently commissioned multi-object spectrograph AAOmega on the 3.9m AAT we have obtained medium- 
resolution near-infrared spectra for 10,500 stars in and around five southern globular clusters. The targets were 47 Tuc, 
M12, M30, M55 and NGC 288. We have measured radial velocities to ±1 km s~^ with the cross correlation method and 
estimated metallicity, effective temperature, surface gravity and rotational velocity for each star by fitting synthetic model 
spectra. An analysis of the velocity maps and velocity dispersion of member stars revealed systemic rotation in four of the 
target clusters. 
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1 Introduction 

Globular clusters are among the oldest objects in our Galaxy 
and their stars provide unique information on Galactic evo- 
lution through cosmic times. These clusters are sensitive in- 
dicators of the Galactic gravitational potential because the 
velocity distribution of the stars in their outskirts is affected 
by the ambient field of the Milky Way. They are believed 
to have undergone substantial dynamical evolution, which 
is affected by processes responsible for the "evaporation" 
of stars (Meylan & Heggie 1997): tidal interaction with the 
Galaxy and two-body relaxation. In most of the cases the 
cluster moves around the Galactic center along very elon- 
gated orbit. Close to the perigalactic position the cluster suf- 
fers the strongest gravitational interactions which can result 
in tidal shocks and lead to the formation of tidal tails. Re- 
cent investigations have found long tidal tails in two low- 
concentration clusters, Palomar 5 and NGC 5466, thus con- 
firming theoretical predictions (Odenkirchen et al. 2001; 
Grillmair & Johnson 2006). The success behind these re- 
sults was secured by deep photometric surveys, most no- 
tably the SDSS project, which allowed cluster membership 
determination with statistical analysis of colors and mag- 
nitudes of stars. Another possibility is offered by spectro- 
scopic observations of individual stars in order to derive 
metallicity and radial velocity, both efficient marker of clus- 
ter members. 

AAOmega, the newly commissioned multi-object spec- 
trograph on the 3.9 meter Anglo- Australian Telescope pro- 
vides an excellent opportunity to measure up to 350-360 ra- 
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dial velocities with a single exposure. Due to its large field 
of view (2 degrees) it can observe stars which are located far 
from the cluster in the plane of the sky. With a velocimet- 
ric accuracy of about 1-2 km s^^, one can easily identify 
escaped stars which have the same or very similar radial ve- 
locities as the host cluster. 

In this paper we present an analysis of radial velocities 
of cluster member stars in terms of systemic rotation for five 
southern globular clusters. A detailed investigation of their 
possible tidal tails has been reported by Kiss et al. (2007). 

2 Observations and data reduction 

The targets were selected from the globular cluster cata- 
logue of Harris (1996) based on the following criteria: we 
preferred nearby, unreddened clusters with large radial ve- 
locities in respect to the Galactic field containing mostly 
disk stars. The latter is very useful to distinguish cluster 
member stars. We also considered recent results on tidal 
tails (in case of NGC 288; Leon et al. 2000), interesting 
dynamical history (Ml 2; de Marchi et al. 2006) or internal 
structures (47 Tuc; Meylan & Mayor 1986) to choose tar- 
gets. We ended up with the following clusters: 47 Tuc, NGC 
288, M12, M30, and M55. By chance we also recorded 
radial velocities of a few stars in the extragalactic globu- 
lar cluster NGC 121 that belongs to the Small Magellanic 
Cloud. 

Our observations were carried out on 7 nights in August, 
2006. Throughout the run the seeing stayed between 1 .2 and 
2.2 arcsec. In total we took spectra for more than 10,500 
stars using the D1700 grating, recording near-infrared spec- 
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Table 1 Total number of stars identified as members, cluster 
radial velocity from Harris (1996) and the measured mean velocity. 



Cluster 


No. 


Vr^d (H96) 


(frad) 




stars 


km s^^ 


km s^^ 


47 Tuc (NGC 104) 


911 


-18.7 


-16.1 


M12 (NGC 6218) 


158 


-42.2 


-40.4 


M30 (NGC 7099) 


129 


-181.9 


-178.1 


M55 (NGC 6809) 


433 


174.8 


171.7 


NGC 288 


123 


-46.6 


-43.5 



tra centered on the Ca II triplet lines. The spectra ranged 
from 8350 A to 8790 A with a resolving power A/ A A = 
10, 500. For each cluster we observed several fields cen- 
tered on the cluster, containing 300-350 stars per configu- 
ration. To reach a signal-to-noise ratio between 50 and 250 
we exposed 60 to 90 minutes in total. For sky background 
measurements and guiding we used 30-40 fibers per config- 
uration. 

The target stars were selected from the 2MASS point 
source catalogue (Skrutskie et al., 2006). We fitted a line to 
the Red Giant Branch (RGB) of each cluster and selected 
stars which matched its color and magnitude. For Ml 2 and 
47 Tuc we selected the lower part of the RGB. Due to the 
sparse field star population around M30 and NGC 288, we 
did not filtered their stars, but used the whole field. The full 
magnitude range of the target stars in X-band was 7 mag 
(from 8 mag to 15 mag) but for a single configuration field 
we limited the brightness range to 3 mag in order to avoid 
cross-talk between the fibers due to scattered light. 

The spectra were reduced with the standard 2dF data 
processing pipeline (drcontml), which extracts automatica- 
lly the wavelength calibrated spectra. Continuum normal- 
ization was done with the IRAF' task onedspec. continuum. 
After this we cleaned the spectra of the remnants of the 
strongest skylines via linear interpolation of the neighbor- 
ing continuum. 

For this study we derived two parameters for each star 
based on its spectrum: radial velocity and the full equiv- 
alent width of the Ca II triplet lines (T,EW = 1^8498 + 
^8542 + W8662)- The latter was determined by fitting a sum 
of a Lorentzian and a Gaussian to the line profiles (Cole 
et al. 2004). Radial velocities were determined in an itera- 
tive way. An initial velocity estimate was provided by the 
line profile fit, which was used to find the best-fit spectrum 
in the extensive spectrum library of Munari et al. (2005). 
The fitted model spectrum was then cross-correlated with 
the observed one, giving the finally adopted velocity. The 
estimated accuracy is about ±1-2 km s^^. 



IRAF is distributed by the National Optical Astronomy Observato- 
ries, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Fig. 1 Heliocentric radial velocities of the member stars ver- 
sus distance from the cluster center. The horizontal lines show the 
mean cluster velocities. Thick lines indicate the catalogued values 
from Harris (1996). The wide distribution of the velocities around 
the mean values is partly due to the intrinsic velocity distribution 
of cluster stars. 

3 Results 

We identified cluster members using the radial velocities, 
equivalent widths and K magnitudes. The selection was ba- 
sed on the Vrad vs. TiEW and YiEW vs. K magnitude cor- 
relations and led to the identifications of 123-911 stars as 
members in the five clusters (see Table 1). Note that these 
sets are more constrained than those of Kiss et al. (2007), 
who only used the radial velocity to determine cluster mem- 
bership. 

First we plot heliocentric radial velocities of the member 
stars as function of distance from the cluster center in Fig. 1 . 
Here the horizontal lines represent the cluster mean radial 
velocities, while on the right-hand side of the plot thick lines 
mark the catalogued values from Harris (1996). While there 
are differences up to 3 km s^^ in the measured means and 
the catalogued values, our samples are more extensive than 
any previous ones for all clusters except 47 Tuc, which may 
explain most of the differences. Because of the excellent 
statistics. Fig. 1 is clearly dominated by the characteristic 
wedge-shaped velocity distribution that can be used to de- 
rive the velocity dispersion profile, which is a sensitive in- 
dicator of the dynamical properties. 

In Fig. 2 we plot celestial positions of member stars for 
47 Tuc and M55, color-coded by their radial velocities. It 
is very prominent that there is a well-defined asymmetry in 
both clusters, which can be interpreted as result of rotation. 

To find quantitative parameters of the systemic rotation 
in these two clusters, we plotted the mean radial velocities 
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Fig. 2 Left panel: stellar positions in 47 Tuc with radial velocities encoded in colors. The circle shows the tidal radius (rt = 42.9')- 
Right panel: the same for M55 (rt = 16.3')- Note the smaller field of view. 
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Fig. 3 Mean radial velocities versus polar angle in 24 azimuthal bins with equal width. We also plotted the best-fit sine curves. 



Table 2 Parameters of the detected systemic rotations: A is the 
peak-to-peak velocity amplitude, while the position angle of the 
rotational axis is given in both equatorial and galactic coordinates. 



Cluster 


A 


PA (eq.) 


PA (gal.) 




km s^^ 


O 


O 


47 Tuc 


6.6±0.07 


130/310±5 


122/302±5 


M12 








M30 


1.5±0.3 


16/196±15 


90/270±15 


M55 


2.4±0.1 


170/350±7 


61/241±7 


NGC 288 


1.4±0.15 


97/277±10 


52/232±10 



versus position angle in Fig. 3. The effects of the internal 
velocity dispersion were minimized by binning the veloci- 
ties in 15 degree- wide sectors of the position angles (a=0° 
corresponding to North, 90° to East). The mean cluster ve- 
locities listed in Table 1 were subtracted for clarity. The si- 
nusoidal shape of the distributions is well recognizable in 
both cases, though it is much clearer for 47 Tuc. For the 
other three clusters similar plots did not reveal obvious pat- 
terns. 

Another useful method to detect azimuthal dependence 
of the radial velocities is stepping an imaginary axis through 



the globular cluster in small angular increments (e.g. one 
degree) and calculate the differences of the mean or median 
radial velocities on both sides of the axis (e.g. Cote et al. 
1995). The resulting curves are plotted in Fig. 4. Based on 
this we can put an upper limit to the rotational velocity pro- 
jected into the line of sight for each cluster by considering 
the amplitude of the curve. For 47 Tuc, the rotational veloc- 
ity is 6.6 km s^^, being in perfect agreement with the results 
of Meylan & Mayor (1986) and Anderson & King (2003). 
The rotational velocity estimates for the other three clusters 
are listed in the second column of Table 2. For M12, the 
curve is clearly non-sinusoidal, so that we did not attempt 
to fit the data. 

To derive the orientations of the projected rotational axis 
in the plane of the sky we fitted a sine wave to the curves 
presented in Fig. 4. The results are listed in Table 2. We also 
converted the equatorial position angles to galactic coordi- 
nates to see if there is a pattern in the directions. Indeed, for 
the four clear detections of systemic rotation we find direc- 
tions being close (within ±20-25 deg) to parallel with the 
Galactic plane, which might be an evidence of non-random 
distribution of rotational axes, but the size of this sample is 
too small even for some statictical approaches. In any case. 
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Fig. 4 Mean radial velocities on both sides of an imaginary axis 
stepped through the clusters. Note the different scale on the y axis 
in case of 47 Tuc. 

we confirmed rotation in 47 Tuc and discovered it in M30, 
M55 and NGC 288. For M12 we will have to clean further 
the sample of the Galactic field contamination before a new 
attempt at detecting rotation is considered. 

4 Future work 

Our data represent a unique source of information on the 
studied five globular clusters. We will separate cluster mem- 
bers from the Galactic field stars by detecting dumpiness 
in the full parameter space of the physical parameters, in- 
cluding radial velocity, metallicity, effective temperature, 
surface gravity. Global kinematics, mass-to-light ratios and 
star formation histories will be constrained from the cleaned 
samples of member stars. We will be particularly interested 
in the age-metallicity relation because a spread in metallic- 
ity can be a sign of prolonged star formation over 2-4 Gyrs 
(Stanford et al. 2006). 

We are specifically interested in investigating mecha- 
nisms that affect velocity distributions in globular clusters 
and, in particular, the tidal tails. Theories to be tested in- 
clude tidal heating of the evaporated stars by the external 
gravitational field (Drukier et al. 1998), the presence of a 
dark matter halo around the clusters (Carraro & Lia 2000), 
and a breakdown of the Newtonian dynamics in the weak- 
acceleration regime (Scarpa et al. 2007). The latter hypoth- 
esis is particularly interesting because modified Newtonian 
dynamics, valid for accelerations below gq 1.2 x 10"^ 
cm s~^, may offer an alternative to the dark matter, with 
far-reaching implications for cosmology. 



For each cluster we aim to measure the velocity disper- 
sion profile of stars, which shows a characteristic decrease 
and flattening at distances from the cluster centre (Scarpa 
et al. 2007). Globular clusters are relatively simple stellar 
systems, with each star's motion dictated by the Newtonian 
gravitational pull of the other cluster members. However, 
any departure from the classical laws of dynamics (such as 
Modified Newtonian Dynamics: MOND, originally intro- 
duced by Milgrom 1983) predicts that the flattening should 
occur at the same absolute acceleration regardless of the 
Galactic environment of the cluster. Hence, we will char- 
acterize the dynamical properties of the clusters as a fun- 
damental probe of Newtonian gravity. Since rotation of the 
cluster can introduce systematic errors in the interpretation, 
we will subtract a smoothed velocity field from the data. 
The residuals will show whether the velocity dispersion de- 
creases at large radii and whether it reverses to an increase 
outside the tidal radius. To disentangle possible breakdown 
of the Newtonian dynamics and conventional dynamic ef- 
fects such as tidal heating, one has to observe clusters of 
very different parameters, including different Galactocen- 
tric distances and determine the actual acceleration at which 
the dispersion profile flattens. This presented sample is a 
good starting point in this direction and while much effort 
will be devoted to model the cluster dynamics, we also plan 
to extend the observational data base with further clusters. 
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